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Abstract 

Two-dimensional heavy-quark QCD is studied in the light-cone coordinates with (anti-) peri- 
odic field boundary conditions. We carry out the light-cone quantization of this gauge invariant 
model. To examine the role of the zero modes of the gauge degrees of freedom, we consider 
the quantization procedure in the zero mode and the nonzero mode sectors separately. In both 
sectors, we obtain the physical variables and their canonical (anti-) commutation relations. The 
physical Hamiltonian is constructed via a step-by-step elimination of the unphysical degrees of 
freedom. It is shown that the zero modes play a crucial role in the self-interaction potential 
of both the heavy-quarks and gluons, and in the interaction potential between them. It is also 
shown that the Faddeev-Popov determinant depends on the zero modes of the gauge degrees 
of freedom. Therefore, one needs to introduce the Faddeev-Popov ghosts in their own nonzero 
mode sector. 



1 Introduction 



One of the promising approaches to problems in QCD is the light-cone quantization [|], [|| . 
The light-cone quantization has turned out to be a useful tool for the perturbative treatment of field 
theories ||, fl. In its extension [9-13] to the nonperturbative domain, one has come to realize that 
careful attention must be paid to the nontrivial vacuum structure of the light-cone quantum field 
theory. Some mathematical aspects of the question with regard to the existence of such vacuum 
states were considered in flTo| , |TT|| . For instance, the light-cone vacuum in the massless Schwinger 
model can be only understood by careful study of the zero modes of the constraints imposed by 
the light-cone frame [12, [ll|. Indeed, it has been conjectured that the dynamics of the zero modes 
in QCD in the light-cone quantization provides the mechanism for the confinement [jl|, g]. 

In the present paper, we continue the quantization of heavy-fermion gauge theory started in our 
previous article [14]. We apply the Faddeev-Jackiw quantization technique to the two-dimensional 
heavy-fermion QED and two-dimensional heavy-quark QCD. Although in this case the quantization 
of the full model (not only in heavy mass limit but also in light mass limit) can be performed, we 
restrict ourselves to heavy- fermions. We carry out the quantization of these models in a light-cone 
domain restricted in its "spatial" directions. It is well known that in such a restricted region one 



has problems with the zero modes [15, [L3fi which, as it was mentioned above, turned out to be the 



most important variables in this case. The quantization of QED and QCD on the finite dimension 
manifolds (circle, torus) were considered in [16, 17, The role of the zero modes in QED using 
the Lagrange approach was studied in fpTsf] . The dynamics of zero modes in the two-dimensional 
QCD employing the light-cone variables was considered in jgg, [21], ^, ]23[]. 



In studying the quantization of the gauge field theories, one is confronted by first-class con- 
straints, and, for this reason, the corresponding gauge conditions should be imposed. A consistent 
canonical quantization formalism for such problems was proposed by Dirac 123] and Bergmann 



[25], and its generalization to fermionic (Grassmann-odd) constraints by Casalbuoni pq| . There is 
another approach to the quantization of the gauge theories proposed by Faddeev and Jackiw [27]. 

The specific gauge theories (QED and QCD) we address in this paper are in terms of the light- 
cone variables where, as we will see, the quantization faces with some constraints involving the zero 
mode variables which require special attention. To examine explicitly the role of the zero modes, we 
consider the quantization procedure in the separated nonzero mode and zero mode sectors. In such 
sectors, we choose special gauge conditions in order to gauge out the nonzero modes and obtain the 
physical variables and their canonical (anti-) commutation relations. The physical Hamiltonian is 
constructed by systematic elimination of the nondynamical variables. 

The paper is organized as follows. In Section 2, the Faddeev-Jackiw technique in terms of 
the light-cone coordinates is applied to the two-dimensional heavy-fermion QED. Despite that the 
results we obtain here are trivial, this example is useful for better understanding what is going on 
in QCD. The Faddeev-Jackiw quantization algorithm is equivalent to the Dirac one, of course, but 
it is sometimes simpler to be employed, especially when the constraints are linear. In Section 3, 
we consider the light-cone quantization of the two-dimensional heavy-quark QCD, where the zero 
modes and nonzero modes of the gauge degrees of freedom are taken into account. It is shown 

(P) a(P) 

that the physical degrees of freedom are the zero modes of gauge fields (A_ in QED and ^4„ in 
QCD), their conjugate momenta and the fermionic variables. We found that the careful elimination 
of unphysical gauge degrees of freedom leads to additional terms in the physical Hamiltonian. Such 
a Hamiltonian is constructed, and, in particular, the interaction potential between heavy-quarks, 
as well as the interaction Hamiltonian between heavy-quarks and gluons are found. It is shown 
that in the case of QCD one needs to introduce the nonzero modes of the Faddeev-Popov ghosts. 
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2 2D Heavy Fermion QED 



In this Section, we are going to consider a simple example which illustrates the Faddeev- 
Jackiw quantization technique [27] for the case of the 2D heavy- fermion QED when one needs to 
take the zero modes of the gauge degrees of freedom into consideration. This will be the first step 
toward the attack of the QCD model, so we can kill two birds by one stone. 

Following the Isgur and Wise [28], the Lagrange density of the 2D heavy- fermion QED has the 
following form 

£ = i^X/KJ^D^ - M\P# + ~F$ 3 (2.1) 

where the Minkowski metric is: diag gn V = (1,-1), 

D„ = 3 M + ieA^ (2.2) 

is the covariant derivative, 

jf =1 . U = r f U(i (2.3) 

and is a given 2- velocity of the heavy- fermions subject to the condition U 2 = 1. The field 
strength tensor is 



03 



d A 3 - d 3 A (2.4) 



and we use the system of units where fi = c = 1. The heavy- fermion limit means that the quantity 
MU^ is greater than any other momenta in the problem under consideration. 

The light-cone coordinates in the two-dimensional space are = (x + ,x~), where 

x ± = -^=(x° ± x 3 ) (2.5) 

The variable x + plays the role of the "time" variable, and x~ is the spatial one. In terms of such 
coordinates, the Lagrange density C becomes 



^ 1 

2' 

e^JjI(U + A^ + U^A + )^ (2.6) 



c = wyr(u + d- + U-d+)^ - + -Fl 



where 



F + _ = d + A_-d_A + , A ± = ±(A ±A S ), 
fe = ^ ■ V± = ±(U ±U 3 ) (2.7) 



The Lagrange density ( |2.1| ) is gauge invariant. This means that the classical theory contains 
"first-class" constraints, and one needs a quantization prescription for systems with constraints such 



as that provided, for example, by Dirac Bergmann |25| , Casalbuoni [26] or Faddeev-Jackiw 

The infrared problems can be regularized by considering the system to be contained in a fi- 
nite volume. We consider, for this reason, the quantization of the theory in the restricted region 
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—L < x < L, and impose periodic boundary conditions for the bosonic variables and antiperiodic 
boundary conditions for the fermionic ones 



^>(x + ,x~ 



-L)=A^(x + ,x +L), 
L) = -<J(x + , x~ + L), \I>(x H 



L) = -$>(x + ,x~ +L) 



[2.1 



We choose the antiperiodic boundary conditions for the fermions in order to avoid their own zero 
mode problem, and treat them as Grassmann-odd variables. The periodic boundary conditions 
for the gauge variables are useful because when integrating by parts the boundary term vanishes. 
Although we choose antiperiodic boundary conditions for the fermions, this does not have to confuse 
anyone because the fermions always appear in a bilinear combination which is a periodic function. 
Each boundary condition implies that there is an additional constraint to be satisfied. In this 
paper, the periodic constraints are imposed explicitly. In [29] the boundary conditions were treated 
as additional constraints. 

To start with quantization procedure, we define the momenta H^, and H^- conjugate to the 
fermionic variables ^ and and the momenta Il-|- conjugate to the gauge degrees of freedom A± 



n± 



~d¥ 
dC 
dA+ 



iwpu-, n 



* — — 







(2.9) 



where the label "r" denotes the right derivative, and the "dot" always means the derivative with 
respect to the "time" x + . 

The zero modes and nonzero modes for the gauge degrees of freedom are defined as 



A 



(P) 



P*A ± , 4« 



Q*A± 



(2.10) 



where P and Q are the projection operators onto the zero mode sector (P sector) and nonzero 
mode sector (Q sector), respectively [13 



P(x,y) 
(P*/)(x) 



1 

2L 
1 

2L 



7TT ' Q(x,y) = S(x -y )-P(x,y) 



f(x)dx 



(2.11) 



The momenta 



n 



(P,Q) - 
± = 



n 



(P,Q) 
i{P,Q) 



conjugate to the variables A^'^ are 



n 



p* 



dC 

dA± 



One then obtains 



n 



(Q) 



Q 



dC 
dA± 



(2.12) 



n 



(P,Q) 



o , n 



(p) 



A 



(P) 



n 



(Q) 



A (Q) 



d-A 



(Q) 



(2.13) 
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The velocities A^'^ can be expressed through the momenta n^ P ' < ^, whereas A^'^ can not. 
Consequently, one has four primary constraints $ = 

f = nf & 
* = < = - wyiu- (2.i4) 
[ <% = % 

In terms of the variables A^f' , n^f' , 11^,, and ^, the Lagrange density C on the constraints 
surface ( p. 14 ) can be rewritten in the Hamilton form 

c = u^A^+il^A^ + u^-h, 

H = Hf + 'Hemi 

h f = -j^-^d-^ -i— n*$fa , 

Hem = -ieU^UZ 1 ^ [U-A { + ] + U^Af ] + E/+aL P) + U+A^} 

+ i ( n ^) 2 + 1 (nL Q) ) 2 + nL Q) a_4 Q) (2.15) 

We neglect the term n^n^ in the Hamilton density Hem because it does not give a contribution 
to the corresponding Hamiltonian. 

We now start the quantization procedure following J27| |. (We do not present here the details of 
this algorithm. The readers can find them in the original paper (2?]] ). To start with, we should 
define the initial set of variables and the corresponding "generalized" momenta a^(C^) which, 
in the case under consideration, are found to be 

Cf = {a^ , nL p) , 4 P) , a {q) , nL Q) , 4 Q) , * , n*} , 

af = {cf , 0,0, Cf, 0,0, C f ,0} , j = l,...,8 (2.16) 



One of the most important objects in the Faddeev-Jackiw approach is the simplectic supermatrix 

f (o; 



fjj, which is defined in general by [| 



Ft „(°) a (0) 

(0) diOj _ ( )eae( dtay 



where e a (ec) ^ s ^ ne Grassmann parity of a(£) 



_ | , if C is Grassmann — even . . 

^ I 1 , if £ is Grassmann — odd 



The simplectic supermatrix corresponding to the set of the variables ( 2.16 ) is block-diagonal 
block diag4 0) (x-,y-) = (A,0,A6(x- -y-),0,Ai5(x- -y~)) (2.19) 

where 
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The appropriate supermatrix fj^ has two eigenvectors with zero eigenvalue Vj = 0) 

y( p ) = {o, o, c ,...o} , 

V^ = {0, ... ,c(x-), 0, 0} (2.21) 

where c is a constant, and c(x~) is a function. This leads to other two constraints 

(ILj,^) {p) = , (2.22) 
S_nL Q) + ie (IL^) (Q) = (2.23) 

After the first reduction procedure, the Lagrange density becomes 

-H F + ie (n**p u + uz 1 a^ q) - \ (n^) 2 - \ (nL Q) ) 2 (2.24) 

where A is a Lagrange multiplier for the constraint ( |2,23| ). 

The constraint ( |2 .22 ) is of the first class (in Dirac classification). Therefore one needs a corre- 
sponding gauge condition. Instead of this, the constraint ( |2.22| ) is to be satisfied on the physical 
state vector. As a result, one can consider the fermionic degrees of freedom, LT^ and as inde- 
pendent ones. Using the constraints ( |2.23| ), one can express, in this case, the momentum in 
terms of the fermionic variables 

n [Q) = -iedZ 1 (IL^) (Q) (2.25) 

where dZ 1 is the operator inverse to 9_, and whose matrix elements in the coordinate representation 
(in the Q sector) are 

G «)(*-- r) = <2^-:£__?C (2.26) 

(One has to keep in mind that the operator dZ 1 is defined only on the Q sector space. Therefore 
( |2.25| ) represents the only solution to the constraint equation ( |2.23| ).) 
The reduced set of variables now is 

Cj 1} = {^ (p) ,nL p) ,^ ) ,nL Q) ,A,^,n,}, J = i,...,7, 

: {C« , , d 1} , , + ieicPcP)™ , C{ 1} , } (2.27) 



4 1} 



V with zero eigenvalue 



The corresponding simplectic supermatrix f-j ' is still singular, but now it has only one eigenvector 



V = {0 , , -d-f(x~) , , f(x~) , -ietf/(aT) , ien*/(x")} (2.28) 

where f{x~) is a function of x~ . This vector does not give any new constraints. Consequently, 
according to p7[| , one needs a gauge condition. 

In the Faddeev-Jackiw quantization procedure, there is only one restriction on how one chooses 
gauge conditions: after imposing gauge conditions, the appropriate simplectic supermatrix should 
be nonsingular. On the other hand, gauge conditions should eliminate the gauge freedom. 
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Let us consider one of the possibilities of how one can choose the gauge condition in the problem 

(2.29) 



considered. Using the constraint (2.23), the term 



(Q) a(Q) 



in the Lagrange density ( fOp can be rewritten as 

A natural gauge condition is chosen to be (see [fHp) 

n G = d^A iQ) = 



(2.30) 



(2.31) 



From ( 2.31| ), it follows that the nonzero mode A^ is an arbitrary function of x + . On the other 
hand, the nonzero mode cannot depend only on x + , otherwise it is the zero mode by definition. 
Therefore, the only solution to the eq.( 2.31 ) is 

A { S ] = (note that A_ / 0) 



(2.32) 



According to ( [2.31 ) and the boundary conditions we have chosen, the term ( 2.30| ) does not give a 
contribution to the Hamiltonian, and can be neglected in the Hamilton density. 
The Lagrange density after the second reduction is 

c m = nL p UL p) +nL Q) iL Q) + n^ + A(9_nL Q) + ie(n^)^) 

+$d„A {Q) - , 



= I(nL p) ) 2 + -(ir_ Q) ) 2 + WF 



(2.33) 



where (3 is the Lagrange multiplier for the gauge condition ( 2.31 ). Now the sets of variables are 

j.i" .ii'' .ir-' . \. ./,• i 8, 

- {Cf } , , Cf , , cL C f + ie{cH\?H Q) , , Cf ) , 0} (2.34) 

and the corresponding simplectic supermatrix has the following block-diagonal form 



a(2) 

,( 2 ) 



f (2) 
Jjk 



A O 

O t B(x-)5(x- 



(2.35) 



where 



/ 





-i 





a_ 








\ 




1 





d_ 


















-d_ 










-ie^> 






-a_ 



























— zell^ 








1 




\ 













1 





/ 



(2.36) 



the matrix A was defined in ( [2.20 ), the matrix O is a zero rectangular 2x6 matrix, and the symbol 
T stands for transposition. 
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(2) 

The simplectic supermatrix f) k is nonsingular, and can be inverted using the Berezin algorithm 
[31] for the inversion of a supermatrix. The result is found to be 



o 

B~ 1 {x~)8{x~ — y~) 



(2.37) 



where 



B-\x~) 



/ 











-dz 1 








\ 










-dz 1 





-iedZ 1 ^ 


iedZ 1 ^ 









dz 1 





dz 2 












dZ l 





-dz 2 





-iedZ 2 ^ 


— iedZ 2 Hy 









-ie^dZ 1 





-ie^dZ 2 





1 




V 





iellydZ 1 





iell^dZ 2 


1 





/ 



(2.38) 



The operator d_ 2 is the one whose matrix elements in the coordinate representation (in the Q 
sector) are 



H(Q)( X - -y-) 



4L 



y~) 2 2L 

3 



(2.39) 



The nonsingularity of the supermatrix (2.37) is consistent with the gauge condition (2.31) we have 
chosen. 

The structure of the supermatrix ( p. 37 ) shows that the physical variables are the fermionic ones, 
11^ and plus the zero modes of Il_ and A_, thereby 



U) 



phys 



{^4 



{P) 



n 



(p) 



satisfying the following brackets 



(P) 



i FJ 



{*(*-) , n*(y-)} 



FJ 



i 

2L' 

5{x~ 



(2.40) 



(2.41) 



These brackets coincide with those obtained in |14| using the Dirac method of quantization for the 
systems with "first-class" constraints. 

The quantization procedure consists of the replacement of the variables u; phys by the corre- 
sponding operators 



phys 



phys 



(2.42) 



which obey the following commutation relations 



[A 



(P) jjCP) 



i 

2L' 



[*(x ), fl*(y )]+ = i6(x -y ) 



(2.43) 



The Hamilton density corresponding to the physical Hamiltonian can be found by solving the 
constraints fl2,23|) , which is equivalent to the substitution of (|2.25| ) into Q2.33|) . This gives 



W phys = I (ttCP)) 2 + f! (g-l * yyru_^y + H F 



(2.44) 
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It should be mentioned that, in the 2D QED with zero modes taken into account, there is no 
real interaction between the gauge field and heavy- fermions. 



Let us discuss the constraint (2.22). As it was mentioned above, it is impossible (at least we do 
not know how to do this) to fix the gauge corresponding to this constraint. We will consider, for 
this reason, the constraint ( |2.22| ) as a strong one, meaning that this constraint should be satisfied 
on the physical state vectors |phys) 

: (IL^) (P) : |phys) = (2.45) 

where : ... : stands for the normal ordering operator. The condition ( [2.45 ) states that the physical 
state vector is chargeless. 

3 2D Heavy Quark QCD 

3.1 Quantization of Heavy Quark QCD 

Here we will consider the generalization of the results obtained in the previous Section to 
the case of non-Abelian gauge theory, in particular, to the two-dimensional heavy-quark QCD. 
We start with the Lagrange density of 2D heavy-quark QCD 

C = %Y^fVU^^f -M^y^f^ff + \fI_FI_ (3.1) 

where the index / stands for the flavor of the quark, f is a color multiplet of quarks with a given 
flavor /, and it is assumed that the quark's mass is flavor independent, 

T> lt = d lt + igA<p a (3.2) 

is the covariant derivative, T a (a = 1, . . . , N% — 1) are the generators of Lie algebra corresponding 
to the group SU(N C ), obeying the following commutation relations 

= if abc T c (3.3) 



rj-id rpb 



with the structure constants f abc being antisymmetric in all indices, and F^_ is the field strength 
tensor 

Fl_ = d+A a _ -d.A a + + g(A + x A_) a , (3.4) 
A% = ±={A%±A%) (3.5) 

Here, for any two "isotopic" vectors, say, T a and Q b , the cross product is defined as 

[f x g) a = f abc r b g c (3.6) 

For simplicity, in what follows, we will consider flavorless quarks. 
The canonical momenta to the variables under consideration are 

n* = iHlfrU- , % = , 

n a(P) = d£_ = A a {P) {P) 
dA _ ~ y ^ > 

U a ^ =Q*J¥i = A a{Q) - d^A a f ] +g(A + x A.) a ^ , 
dA_ + 

n a(P,Q) = Q ^ ^ 
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The velocities A a ^ P ^ can be expressed through the momenta n"^' , and one has the following 
primary constraints $ = 

( ^ a ( p -o) = n a(p ' Q) 
* = < = - i^yju- (3.8) 
{ <% = % 

On the surface of these constraints, the Lagrange density ( |3.1| ) in the Hamilton form reads as 



where 



c = i^ p) A a } p) + if} Q) A a } Q) + n*^ - n , (3.9) 



W = ~Hf + ^Gj 



Hg = \ (n a j p)2 + n a j Q)2 ) + Yf} Q) d^ Q) - g n a ^ (A + x a_ 



va(P) 



-igU^T a ^> (A a } P) + (3.10) 

Using the similar arguments from the previous Section, we neglect the term n_ P ^n_^ in the 
Hamilton density TLq- 

The sets of the initial variables an d the "generalized" momenta a^°'(C) are 

C (0) = | A a(P) n a(P) A aiP) A aiQ) ^ ^(Q) ^Q) ^ ^ ^ ^ j 



x (0) 



={d° } , 0, 0, Cf , 0, 0, C^ 0) , 0} , j = l,...,8 (3.11) 



The corresponding simplectic supermatrix has two eigenvectors with zero eigenvalues 



V^ = {0, 0, c a ,...0} , 

= {0 , ... , c a {x-) ,0,0} (3.12) 

leading to two other sets of constraints 

ig (n*T a ^) (p) + g (A- x n_) a(p) = , (3.13) 

d _jfiQ) + g ( A ( p ) x n W) ) a + g (a (q) x n_) a + i g n*T a y = o (3.14) 



Equation ( 3.13| ) tells that the total classical color charge in the system is zero. 

The Hamilton density TLq on the constraint surface (|3.8|), ( |3.13| ), and ( |3.14j ) becomes 

n G = \ (n a ^) 2 + 1 (n a ^f - i g ^u^(A a ^ + a«®) (3.15) 

If one adds the constraints ( [3.13 ) and ( p. 14 ) with corresponding Lagrange multipliers to the La- 
grange density ( |3.9j ), this does not give any new constraints. Therefore, one needs gauge conditions. 
Let us choose them, like in QED, in the form 

n a G = d-A a _ {Q) = (3.16) 
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Using the same arguments as in the previous Section, one obtains 

A <Q) = 0> ( no te that At + 0) (3.17) 
The Hamilton density Tic then can be sufficiently simplified 



It can be shown, using the constraints ( 3.13 ), that the last term in the above expression does 
not give a contribution to the Hamiltonian Hq. Keeping this in mind, one can neglect it in the 
Hamilton density. Therefore, one can write 

Ha = 1 -(u a ^)\ 1 -(n a ^) 2 (3.18) 

where the variables have to be eliminated using the Gauss law constraints 

W GL = d„Tf} Q) + g (A^ x H W) ) a - g (a^ x nL P) ) a + igp a = (3.19) 

where 

p a = (n*T°^) (Q) (3.20) 
In deriving fl3.19| ), we have used the relation 

(A_ x n_) a ( p ) = {A {P) x U (P) ) a 

which holds on the surface of the gauge condition ( |3.16 ) (or fl3.17|) ). 

f.(?" 



(2) 

Now the corresponding simplectic supermatrix ' is nonsingular, and has a similar form as in 



2D QED with 



{A« p \^}fj = ±6* 



{n»(x") , *(y-)} F j = S(x--y-) (3.21) 
Therefore, one can consider the set of the variables 

;P hys = |(^( p ) ; n " (P) , II*, (3.22) 



lO 1 



as physical ones. The quantization procedure consists of the replacement of the variables w phys by 
the corresponding operators 

w Phys ^ ^Phys ( 3 23 ) 

which obey the following commutation and anticommutation relations 

[A a _ (P \ frl (P) ]_ = ^5 ab , (3.24) 
[*(x"), n*(2T)]+ = iS(x~ - y~) (3.25) 
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3.2 Physical Hamiltonian 

To obtain the physical Hamiltonian, one has to solve the constraints ( |3.19| ) in order to 
express the momenta in terms of physical variables. Although the eq.( 3.19| ) is the one for the 

nonzero mode of the momenta IT^f , its left-hand-side contains a term which is a zero mode (the 
product of the zero modes does not have the nonzero mode). Thus, one has to be careful when 
such an equation has to be solved. 
Let us define the function 



z a { x ~) = v ab (x-)n aw \x 



a(Q), 



where 



(3.26) 
(3.27) 



V ab (x-) = (exp( 5 i?x-)) afe , (V 6 = f acb A c S P) 
is a unitary matrix. The Hamilton density ( p. 18 ) can be expressed through the function Z a (x~) 

H G = ^(n a _ (P) ) 2 + Uz a (x-)) 2 (3.28) 



From ( |3.19D , the equation that the function Z° 

d-Z a (x-) = (V{x~)L{x-)) a , 

L a {x~) = g (a (P) x U (P) Y - igp a {x- 

The solution to the eq. fl3.29 ) has the form 

Z a ( x ~) = (VL) a{p) x- + dZ 1 (VL) a{Q) (x-) 



(3.29) 
(3.30) 

(3.31) 



We wish to stress out the existence of the linear term in ( |3.31| ). As we will see later, this term 
will give a contribution to the physical Hamiltonian where the momenta are involved. 

To obtain the zero modes and nonzero modes of the expressions involved in (3.31), we remark 



that 



L a(P) 



L a(Q) 



(3.32) 



Substituting these expressions into (3.26), one obtains the function Z a (x ) in terms of the physical 
variables uj^ 8 defined by (gjg ) 



(3.33) 
(3.34) 



Z\x~) =g{v ac ^ (A iP) x rt P) y-iK a ^}x- 

-iK b ^-iV bc ^p c (x-)} 



+ J r 1 {(r(«'( I -)(/ ) xn (p) 

where 

K a (x-) = V ab{ - Q) {x-)p b {x-) 
The eq.( |3.33| ) makes possible to construct the physical Hamiltonian H phy: 



^-phys 



n phys dx~ , 

—L 

H F + \{lf^?+ l -{Z«{x-)f 



(3.35) 
(3.36) 
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where Z a (x ) is given by ( 3.33| ) and can be presented in the following form 



Z a {x~ 

B a(P 

G a (x~ 

r a{Q) (x- 



G a(P)_ iK a(P) 

G a{ - Q \x-)-iV a{ - Q) {x-) 
V ab (x-)(Ai P) xU^) b = 
K< Q \x-) + V ab ^p b 



(3.37) 



V ab (x-)R bc U c S P) 



(3.38) 



Establishing the necessary tools, we can now consider the interaction Hamiltonian. Due to 
eqs.( |3.36| ) and ( |3.38D , the interaction Hamiltonian corresponds to the last term in ( 3.36| ) 

Q 2 f L 

H int = y J_ L H mt (x)dx (3.39) 

where 

H- mt {x) = B a{p) B a{p) x 2 + 2xB^ p) d- 1 B a{ - Q) {x) + d~ 1 B< Q \x)d' 1 B< Q \x) (3.40) 
One can present the interaction Hamiltonian in the following form 

H mt = H%°° + + (3-41) 



where Hf^° is the part of interaction Hamiltonian which describes the self-interaction of the gauge 
zero modes in the sector of the pure Yang- Mills fields 



rrzero 
-^int 



9 



—G a ^G a ^+ [ L d~ 1 G a ^\x)d- 1 G< Q \x)dx 

3 J-L 



-2G a ^ 



L ,x 2 L 2 



)G< Q \x)dx 



(3.42) 



-L 2 6 

the term //'^ ear is the part of interaction Hamiltonian which is linear over the quark color charge 

p a 

~2L 3 rL 



linear 
int 



-w 



^ G a(P) K a(P) + f d -^ G aiQ) {x)d -l V a{Q) {x)d 
3 J-L 

, (G a{ - p ^V a{ - Q) {x)+K a{ - p ^G a{ - Q \x)) dx 



L x 2 I? 



( 9 
L 2 



(3.43) 



and the term describes the quark-quark interaction 



H, 



int. 
quad 



int 



2L 3 



_ 2K a{P) 



K a(P) K a{P) 



+ 



d- 1 v a{ - Q \x)d- 1 v a{ - Q \x)dx 



L ,x 2 L 2 



)V a{Q) {x)dx 



-L 2 6 

In deriving these expressions, we used the following properties of the operator <9 _1 
{d- 1 f[ Q) ){x)i Q \x)dx = - ( L {d- x i Q) ){x)f[ Q \x)dx, 



(3.44) 



L 

1 L 2 



d 1 x = 

2 6 



(3.45) 
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For the Abelian case, the only term that survives in (CT) is F^ ad , given by (pOl ) with if a = 
and = (ILi>^) (Q) . This gives the result obtained in the previous Section. 

Let us now return to the constraint (3.13). Using the gauge condition ( 3,16| ) (or ( 3,17| )), it can 



be rewritten as 

i (n^r a ^) (p) + (a {p) x ni p) ) a = o (3.46) 

This constraint is of the first class. We do not know how to fix the gauge corresponding to this 
constraint. Therefore, as in the previous Section, we will consider the constraint ( [3.46 ) as a strong 
one, meaning that it should be satisfied on the physical state vectors |phys) 

: (i (II^T^ ) (p) + (a (P) x IlL P) ) a ) : |phys) = (3.47) 

where : ... : stands for the normal ordering operator. 

3.3 The group SU(2) 

Let us examine the results obtained for the SU(2) group. In this case, the structure constants 
jabc coincide w ith the antisymmetric Levi-Civita tensor e abc 

f abc = e abc ; a = 1)2j3) ( e 123 = l) (3.48) 
The matrix V ab (x~) can be calculated and is found to be 

/ T>2\ab P a ' ) 

V a \x-) = d ab + V l (x-)^- + V 2 (x-)^-, 
Vi(x~) = 1- cos (gx~D) , V 2 (x~) = sm(gx~D) (3.49) 



D = J A a _ {P) A a _ iP) , (R"Y b = A a _ {P) A b S P) - 6 ab D* (3.50) 



where 



From these expressions, one can explicitly find the zero modes and nonzero modes of the matrix 
V ab (x~) involved in the physical Hamiltonian 

I r>2\ab f? a ^ 

V*W(x-) = ui(x-)^-+« a (x-)^- > 

I n2\ab 

yab { P) = 6 a b _ Ui{0) («_) 

D z 

sin( gDL) 

u\(x ) = jjj-^ cos(gDx ) , U2(x )=s'm(gDx ) (3.51) 

Other quantities that are involved in the problem are 

( r>2\ab nab 

dZ^HQ)^-) = Vl{x -)^_ + V2{x -)^_ ^ 

-\ 1 -\ _s'm( gDL) _ 1 _. 

vi(x ) = ^U2{x ) + x — — , v 2 {x ) = — ui(x ), (3.52) 

gD gDL gD 

->2\ab pafe 



d-JV^\x-) = zi(x-)^- + z 2 (x-)^, 
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zi(x 
z 2 (x 

G a(P 

G a(Q) (x 
K a (x 



1 . . Ax") 2 L 2 ^sm(gDL) 
u 1 (x-) + C 



g 2 D 2 
1 

g^D 2 



u 2 (x ) + 



s ' m 9DL abu b(P) 
gDL 



( R ab 
D | m(x~) — - u 2 {x~) 



2 6 ; gDL ' 
x~ s'm(gDL) 
1)D gDL ' 

{R 2 ) ab ' 



(3.53) 



D 2 i - ' 



^ (R 2 ) ab R ab \ 

Mx~)^r- + ^(x-)— J P \x-) 



(3.54) 



The part Hf^° of the interaction Hamiltonian can be simplified and presented in the following 
form 



rrzcro 
^int 



sm(gDL) Y 
gDL J 



a{ P)(R 2 r b 5 ( p) 

L> 2 



This term describes the self-interaction of the gluons in terms of their zero modes. 
The result for the part of the interaction Hamiltonian Hj^ ear is 



rj-lincar 
^int 



-ig 2 D\ I Ai(z)IT 



a(P) (R 2 T b J, 

D 2 



L , D , Dab > 

p b {x)dx+ / ^A 2 (x)n a _ (P) ^-p 6 (x)dx 



(3.55) 



(3.56) 



where 



x sm(gDL) cos(gDL) 
Xi(x) = — — — cos(gDx) H — sm(gDx), 



A 2 (x) 



yD yDL 
x sin(y.DL) 
~g~D~g~DL~ 



g 2 D 2 



1 



sin^Dx) - 3^ cos(gDL) 



sin(gDL) 

gDL 



- cos(gDx) (3.57) 



This term describes the interaction of the heavy-quarks with the zero modes of the gauge degrees 
of freedom. 

The part H^ ad of the interaction Hamiltonian, which describes the self-interaction between 
heavy-quarks, has the following form 

rL r L / 



jj-quad _ 2 
-"int - 9 



L J-L 



A 1 (x,y)p a (x) [ -^p b (y) + A 2 (x,y)p a (x)—p b (y)\dxdy 



W j\ j\ P a (*) (t ab + P b (y)H^\x - y)dxdy 



(3.58) 



where 



Ai(x,y) = ^ (ui(x)ui(y) + u 2 (x)u 2 {y)) 



+ 



2L 



x 2 L 2 

(- — ) cos(gDy) - u 2 (x) sin(gDy)) 



+ [x ^ y] 



+ ^ [(«i(^) cos(gDy) - u 2 (x) sm(gDy)) + (x <-> y)] F (Q) (x - y) 



sin(yDL) 
IgDL 



(cos(gDx) + cos(gDy)) H {Q) (x - y), 



(3.59) 
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M{x,y) 



L 
6" 



{u-i_{x)u2(y) - u 2 (x)u 1 (y)) 



1 

2L 



x 2 L 2 

(— — ) (u 2 (x) cos(gDy) + m(x) sin(gDy)) 



[x <-> y] 



o [(^0*0 cos(gDy) + ui(x) sin(gDy)) - {x y)] F (Q) (x - y) 



+ S1 ^ L L) (sin(yDx) - sin( 5j Dy)) - y) 



(3.60) 



where H^\x — y) is given by ( 2.3S| ). 

Prom ( p6|) and ( p5|) , one can find the total Hamiltonian corresponding to the pure Yang-Mills 
fields 



tt _ rT[ a{P) r fab T Tb(P) 

-Hpure YM — LLL_ Lr 11_ 



(3.61) 



where 



2\ab 



smgDL y\ (R 2 ) 
gDL J J ~D2" 



(3.62) 



The Hamiltonian H pure ym is quadratic in the momenta and sufficiently nonlinear in the 



coordinates A 



a(P) 



3.4 The Faddeev-Popov Determinant 

The constraints ( |3.19|) and gauge conditions ( |3.16|) are the only local constraints in the 
problem. Although these constraints are linear over the variables and A^!P\ their Poisson 

bracket is nontrivial and depends on the physical variables A_ . This might give a nonzero 
contribution to the effective interaction potential. To find such a contribution, we should calculate 
the Faddeev-Popov determinant which appears due to the constraints ( [3.16 ) and (3.19). One should 



introduce the Faddeev-Popov ghosts provided that the Faddeev-Popov determinant is not trivial 
(it may depend on the field variables). On the other hand, the Faddeev-Popov determinant is 
important for the unitarity of the S-matrix. 

In order to find the Faddeev-Popov determinant, one should consider the Poisson bracket be- 
tween the first-class constraints fl3.19| ) and the corresponding gauge conditions ( |3.16 ). The result 
is found to be 

{n a G (x) , n b GL (y)} PB = -d x V a }{x)D®\x-y) , (3.63) 
V^{x) = 5 ab d x + gR ab , R ab = f acb A°S P) (3.64) 

where 

D (Q)( X -y)= S ( X -y)--L (3.65) 

is the (^-function in the Q sector. (For simplicity, we omitted the index "-" in the coordinates x 
and y.) 

The Faddeev-Popov determinant is defined as 

A[A] = Det 1 / 2 ^, 9ff} PB \9=o (3-66) 
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where is the set of first class constraints and the corresponding gauge conditions. In our case, 

(p) 

the Faddeev-Popov determinant depends only on the variables and can be presented as 



A[A^>] = Det 1-^ j = Det ^— j (3.67) 

In order to normalize the Faddeev-Popov determinant, we have introduced the operator in the 
denominator of ( 3.67| ) . Using the well known property of the determinant of a matrix 

DetM = exp (Tr LnM) 

one obtains 



A[A ( _ J ] = exp l^Tr Ln-^-J (3.68) 

Taking the derivative of the both sides of ( 0-68 ) with respect to the coupling constant g, one 
obtains an ordinary differential equation for the Faddeev-Popov determinant 



(TrVZ 1 (g)R) A[A a _ {P) ] = (TrG(\g)R) A[Al (p) ] (3.69) 



dA[A a _ {P) ] 
dg 

where G ab (x,y\g) is the Green function of the operator T>°^- The solution of this equation has the 
following form 

AL4 a _ (P) ] = exp (J* dg'TrG(\g')R^) (3.70) 
The determinant ( 3.67| ) is equivalent to the Faddeev-Popov term in the effective action 

S FP = J dx + dx-^ 1 {x)d-V ah {x)c\x) (3.71) 

where c b and c a are the ghost and antighost fields, respectively. The operator T>°^{x) depends only 

on the zero modes A a ^ P \ Therefore, only the nonzero modes of the ghost and antighost fields give 
contributions to the Faddeev-Popov action, meaning that the Faddeev-Popov ghosts are necessary 
only in their own Q-sector. One should then make the following substitution 

c a _ > c a ( < 9) c a ^ c a (Q) 

and the Faddeev-Popov action becomes 

S FP = J dx + dx'c a ^ Q \x)d-V ah (x)S Q \x) (3.72) 

where and are the nonzero modes of the Faddeev-Popov ghosts. 

Let us compute the Faddeev-Popov determinant. As it follows from (|3.7C| ), one needs to know 
the Green function G ab (x,y\g). According to ( |3.72| ), it satisfies the equation 

d x G ab (x, y\g) + gR ac G cb (x, y\g) = 6 ab D^ (x - y) (3.73) 

The matrix R ab does not depend on x. Consequently, the Green function G ab (x, y\g) is translational 
invariant: G ab (x, y\g) = G ab (x — y\g) = G ab (z\g), and ( |3.7lD can be rewritten as 

d z _G ab {z\g) + gR ac G cb (z\g) = 5 ab D^ Q \z) (3.74) 
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This equation is an inhomogeneous first order differential equation, and its solution can be presented 
in the following form 



G ab (z\g) = U ac (z\g) (v cb (\g)D^) (P) z + U^z^d' 1 (v cb (\g)D^) W> (z 



(Q) 



(3.75) 



where the matrix V ab (z\g) was defined by ( |3.27| ) , and the matrix U (z\g) is its inverse one. Sub- 
stituting (|3~7^ ) into (|3?70|) , one obtains 



A[A a _ (P) ] = exp (2Lj° dg'd- 1 {y ah (\g')D^ Q) {z)\ z=0 R 
The integrand in the right-hand-side of eq. ( 3.76| ) can be presented as 

d- 1 (y ab {\g')D^) iQ) (z) = 

= V ab ^{g')G^\z)+d- 1 (y ab ^\\g')D^) (Q) (z) 



(3.76) 



where the function G^\z) = 4p — ^ was defined by ( |2.26 ), Thus 



A[A a _ {P) ] = exp (2L j°dg'd- 1 (v ab ^(\g')D^ 



(Q) 



R 



ba 



(3.77) 



(3.78) 



This is, so far, the most general result we can obtain for the Faddeev-Popov determinant. We 
cannot go further without specifying the gauge group. Therefore, for simplicity, in what follows 
we will consider the group SU(2). Using the representation ( |3.51 ) for the matrix y ab (<2) and the 
definition of the operator d , one obtains 



A[A 



a(Ph 



exp 



4DL ["dg'd- 1 (u 2 (\g')D^) iQ) (z)\ z=0 



d- 1 (u 2 (\g')D^(z)) iQ) \ z=0 



1 



2g'DL 



sin 



-(g'DL) \ 
g'DL ) 



Thus 



A[A a _ {P) ] = exp 



gDL 



SIM 



dx 



(3.79) 
(3.80) 

(3.81) 



The Faddeev-Popov determinant depends on the zero modes A°^ P ^ through the upper limit 
in the integral involved in ( 3.81 ). This means that one does need to introduce the Faddeev- 
Popov ghosts (at least for the group SU(2)), and there is an additional term to the self-interaction 
Hamiltonian of the gluons. 

The integral that appears in the right-hand-side of ( |3.81 ) can be presented in the following form 



smx 



dx 



sm z 



+ (Ci(z) - 7 - In z) 



where Ci(z) is the cosine integral special function, and 7 is the Euler constant. Therefore 



A[A a } P) ] = exp<{2(l 



H 1 



sm 



(gDL) \ 

gDL J 



2 (Ci(gDL) - 7 - In(gDL)) 



(3.82) 



(3.83) 
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4 Conclusions 



We have considered the light-cone quantization of the two-dimensional heavy-quark QCD, 
explicitly taking into account the zero modes contribution of the gauge degrees of freedom. We have 
imposed the periodic boundary conditions for the gauge degrees of freedom and the antiperiodic 
ones for the fermions. As ordinary QCD, this model is gauge invariant, meaning that there are 
unphysical degrees of freedom in the problem. We have used the Faddeev-Jackiw algorithm to 
quantize the theory . In order to see the role of the zero modes explicitly, we have considered the 
gauge conditions which are needed only in the nonzero mode sector. We obtained the physical 
variables (coordinates and their conjugate momenta) and the corresponding (anti-) commutation 
relations. We found that the physical variables are the zero modes of the "spatial" light-cone 
gauge degrees of freedom A°^ P ^ and n"^, and the fermionic variables \P and II Solving the 
constrains and the gauge conditions in order to eliminate the unphysical gauge degrees of freedom, 
we constructed the physical Hamiltonian. In the elimination of the unphysical variables mentioned 
above, we found that the expression for the momenta n_ contains an additional term that is 
linear in x~ (see ( |3.31| )). Such a term gives a contribution to both the self- interaction and mutual- 
interaction potentials of the gluon fields and quark fields. We found that one needs to introduce 
the Faddeev-Popov ghosts in their own nonzero mode sector. The Faddeev-Popov determinant was 
calculated, and it was found that it depends on the zero modes A°^ P ' , giving a contribution to the 
self- interaction Hamiltonian of the gluons. 
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